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Abstract 

We perform a fully relativistic calculation of the D(e,e'p)n reaction in the impulse approxima- 
tion. We employ the Gross equation to describe the deuteron ground state, and we use the SAID 
parametrization of the full NN scattering amplitude to describe the final state interactions (FSIs). 
We include both on-shell and positive-energy off-shell contributions in our FSI calculation. We 
show results for momentum distributions and angular distributions of the differential cross section, 
as well as for various asymmetries. We identify kinematic regions where various parts of the final 
state interactions are relevant, and discuss the theoretical uncertainties connected with calculations 
at high missing momenta. 
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I. INTRODUCTION 



Exclusive electron scattering from the deuteron target is very interesting by itself, and 
also as a very relevant stepping stone towards understanding exclusive electron scattering 
from heavier nuclei. The -D(e, e'p)n reaction at GeV energies allows us - and requires us - to 
carefully study the reaction mechanism. It is necessary to consider final state interactions 
(FSIs) between the two nucleons in the final state, two-body currents, and isobar contri- 
butions. Of these, the FSIs can be expected to be the most relevant part of the reaction 
mechanisms at the GeV energy and momentum transfers relevant to the study of the tran- 
sition from hadronic to quark-gluon degrees of freedom. For some recent reviews on this 
exciting topic, see e.g. 

naa. 

Even though the deuteron is the simplest nucleus, and it has been the subject of consid- 
erable attention for a long time, there are several open questions: Is it possible to exper- 
imentally determine the high-momentum components of the deuteron wave function even 
though the nuclear wave function or the momentum distribution are not observables? Will 
conventional nuclear physics break down or become too cumbersome at some point, and will 
a description involving quarks become necessary? Are there any six-quark admixtures in 
the deuteron wave function, and do they have an unambiguous experimental signal? What 
influence does short range physics have on conventional wave functions, and how can this 
influence be removed to find an effective potential Viowk, provided one is interested in low 
energy scenarios only P, ^? Relativistic wave functions are available for the deuteron. 
And while the calculational effort is still considerable, it can be managed without resorting 
to super computers. The interest and importance of the D{e,e'p)n process is reflected in 
the fact that recently, a deuteron benchmarking project has been started to investigate the 
differences between various calculations that are based on non-relativistic wave functions 



0. 



Anything that we learn about the reaction mechanism of the D{e,e'p)n reaction has 
implications for heavier targets or experiments where the deuteron is used as a lab. Some 
examples for the latter are the measurement of the neutron magnetic form factor by measur- 
ing a ratio of D{e,e'p) and D{e,e'n) cross sections. This allows for a significant reduction 
in the model dependence of the extracted form factor, but some theoretical input is still 
required IS]. Another example for using the nucleus as a lab is color transparency. While 
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meson production from nuclei 
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12| is the main thrust of color transparency inves 



In 



tigations, color transparency in (e, e'p) reactions is very interesting and topical, too 
order to study color transparency, one first needs to establish a firm understanding of all 
the conventional nuclear effects. 

Several experiments on deuteron targets have been performed in the last few years, both 
at Jefferson Lab and at MIT Bates, and these data have either been published recently or 



are currently under analysis 
experiments at Jefferson Lab 



HQ, 



171]. There are also new proposals for D(e,e'p) 



. Apart from these exciting new data, there are very in- 



leresting open questions posed by the data that have been available for some years now 
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20 . 121I . |22| . Regardless of the momentum transfers involved, there has been a discrep- 



ancy between data and calculations at low missing momentum. We discuss calculations in 
kinematics relevant to the new experiments as well as the low missing momentum puzzle in 
our results section. 

The experimental activity in (e, e'p) reactions on the deuteron and other light nuclei 
has been matched by theoretical efforts. These calculations typically are performed using 
Glauber theory, the generalized eikonal approximation 
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251], or a diagrammatic ap- 



271 ]. Even a second-order correction to the 



proach [26|, although there are rare exceptions 
eikonal approximation has been suggested recently [28|. Many of these calculations focus 
on the differential cross section only, and use just the central part of the NN scatte ring am- 
plitude. Currently, almost all calculations for D{e, e'v)n reactions are unfactorized [29[, but 
factorized approaches are used for heavier targets 30|]. A common feature first introduced 
in Glauber calculations is the assumption that the momentum transfer in the rescattering of 
the two nucleons is purely transverse. This has consequences both for the profile function, 
and the argument for which the NN scattering amplitude is evaluated. 

In this paper, we present a new calculation with several important features: we use 
a fully relativistic formalism, and describe the ground state with a solution of the Gross 
equation [311]; we include all parts of the nucleon-nucleon scattering amplitude, including all 
the spin dependent parts, and use a realistic, modern parametrization; the only approxima- 
tion that we make is to neglect the negative energy states, as discussed at the end of the 
first section. This new calculation can be used at all energy and momentum transfers, pro- 
vided that an appropriate full set of pn scattering data is available. This paper is organized 



as follows: First, we review the theoretical framework for our calculations, then, we show 
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(a) 



(b) 



FIG. 1: Feynman diagrams representing the impulse approximation. 

numerical results for the cross section and asymmetries. We conclude with a summary and 
outlook. 

II. THEORETICAL FRAMEWORK 



The Feynman diagrams representing the impulse approximation are shown in Fig. [H 
Figure [1^ represents the plane wave contribution and Fig. [TJo represents the contribution 
from final state interactions. The plane wave contribution to the current matrix element is 
given by 

{Pisi; P2S2I J^w \P^d) = -uipu si)r^{q)GoiP - P2)rl^{p2, P)u^iP2, S2) , (1) 

where the target deuteron has four-momentum P and spin A^, the final proton has four- 
momentum pi and spin Si and the final neutron has four- momentum p2 and spin 82- The 
single-nucleon propagator is 

Goip) 



J ■ p + m 



rn? — p"^ — IT] 

and the current operator is chosen to be of the free Dirac-plus-Pauli form 
The deuteron vertex function with nucleon 2 on shell can be written as 

.P-i^AP) 



(2) 



(3) 



TX,{P2.P) = 9l{piP2-P)l-ixAP)+92{piP2-P)- 

- (g3{plP2 ■ P)i ■ ^xAP) + 9a{pIp2 ■ P)- 



m 



P-^\diP)\ 1-Pi + m 



m 



m 



C, (4) 



where Pi — P — P2, P — ^{pi — Pi) — ^ — P2-, C \s the charge-conjugation matrix and is 
the deuteron polarization four-vector. The invariant functions Qi are given by 



gi{pl,P2 ■ P) 
92{pI,P2 ■ P) 

g3ipl,P2 ■ P) 

94{pI,P2 ■ P) 
where 



2Ek - M, 



2Ek - Md 



Stt 



1 / 3 TTT- 



— m(A;) H + \ --rvAk) 



3 mEk 

IGtt A; 

2 

m 



87rMd 



Ek + m 



V2k^ 



k 



Vs{k) 



k 



{P-P2? 2 
-P2 



p2 



is the magnitude of the neutron three-momentum in the deuteron rest frame and 

Ek = + m? . 



(5) 
(6) 
(7) 

(8) 
(9) 

(10) 



The functions u{k), w{k), Vs{k) and Vt{k) are the s-wave, d-wave, singlet p-wave and triple 
p-wave radial wave functions of the deuteron in momentum space. The radial wave functions 
are normalized in the absence of energy-dependent kernels such that 

(2^ 



^0 



[u\p) + w\p) + vlip) + vHp)] = 1 



(11) 



For convenience, the spectator deuteron wave function can be defined as 

V'A,,..(P2, P) = G'o(P - P2)rl,(p2, P)U''{P2, S2) . 

We choose to normalize this wave function such that in the deuteron rest frame 

(Pp2 m 



(12) 



S2 



Y E 



-^Xd,S2{P2: P)l ^Xa,S2{P2: P) = 1 



(13) 



P2 



which is correct only in the absence of energy-dependent kernels. The plane wave contribu- 
tion to the current matrix element can then be written as 



{PlSi;P2S2\ Jpw \P^d) = -U{P1, Si)T''{q)^Xd,S2{P2, P) 



(14) 
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The contribution from final state interactions represented by Fig. [T]d requires an integra- 
tion for the loop four-momentum k2 which involves both the deuteron vertex function and 
the pn scattering amplitude. An equivalent approach is to formulate the problem in terms 
of the Spectator or Gross equations 3l| where the equations for the scattering amplitude 
and vertex function are rewritten such that one particle is always taken to be on mass shell. 
This approach is manifestly covariant and has been successfully applied to elastic electron 



scattering from the deuteron [3^. Using this approach, the contribution of the final state 
interaction to the current matrix element is given by 

77— r— Ma(Pl, Si)Ufe(p2, S2)Mab-cd{j>l,V2; h) 

xAj,{k2)Tljh,P), (15) 
where M is the pn scattering amplitude, 



2m 



(16) 



is the positive energy projection operator and the Dirac indices for the various components 
are shown explicitly. Since the single-nucleon propagator can be decomposed as 



Go{p) 



m X ^ 



p 
m 



and 



(fT5|) can be written as 



En 



1 



m(p, s)u{p, s) ^ v{-p, s)v{-p, s) 
_p<^- Ep + ie 

A+{p) A- 



p^ + Ep- ie 
-P) 



Ep + ie p^ + Ep — ie 
V 



p^ - Ep + ie 



-tTrSip" -Ep) + - 



En 



(17) 



(18) 



{PlSl]P2S2\ JpSI \P^d) = {PlSl]P2S2\ Ja I -P^d) + (piSi ; 1 |-PAd) + (piSi ; 1 Jc \P^d) 

(19) 

where 

d?k2 m m 



{piSi;P2S2\J^\PXd) = / 



(27r)3 Ek, Ep+,_ 



k2 



5(pu j^j,_Ek,- Ep+^_k,) 



XUa{pi, Si)Ub{p2, S2)Mah-cd{Pl,P2] k2)Ud{k2, (12) 

X (A+(P + q-fc2)r^(g)^A„.2(^2,P)), , (20) 
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(a) (b) (c) 

FIG. 2: Diagrams representing {piSi;p2S2\ Ja \PXd) , {PiSi;P2S2\ Jb \P^d) and 
{piSi;p2S2\ Jc \PXd) ■ The insertion of a cross on an internal propagator indicates that the 
particle has been placed on the positive-energy mass shell. The insertion 0+ indicates that the 
positive-energy principal value part of the propagator is used and the insertion of O— indicates 
the negative-energy part of the propagator is used. 



XUaiPl, Si)Ub{p2, S2)Mab;cdiPl,P2; k2)Ud{k2, CTa) 

X (A+(P + q-fc2)r^(g)^A„..(A;2,P)), (21) 



and 



I , X - [ d?k2 mm 1 

{PlSuP2S2\J^\PXd) = > / 77— T^ — ^^ — 

XUaiPl, Si)Ub{p2, S2)Mab-cd{pl,P2; k2)Ud{k2, (T2) 

X (A-(fc2 -q- P)T'^{q)^x„^,ik2,P))^ • (22) 

These contributions are represented by the diagrams in Fig. |2J 

Equation ( 120|) . represented by Fig. [2ti, has all four legs of the pn scattering amplitude 
on mass shell. For this case, the scattering amplitude can be parameterized in terms of five 
Fermi invariants as 

Mab;cd = ^sis, t)SacSbd + ^v{s, t)'^ac " Ibd + ^t(s, t)'^ac i(^f^u)bd 

+J'pis, t)^llld + ^a(s, t){^'l)ac ■ {i'i)m (23) 



where s and t are the usual Mandelstam variables. The calculation of the invariant functions 
from helicity amplitudes is described in Appendix |X1 We use the helicity amplitudes available 
from SAID as input for our calculations [33]. For this calculation we have constructed a 
table of the invariant functions in terms of s and the center of momentum angle 6. The 
table is then interpolated to obtain the invariant functions at the values required by the 
integration. 

An alternative two-dimensional representation of the scattering amplitudes is in terms 
of the Saclay amplitudes. In this case, the scattering amplitude as an operator in two- 
dimensional spinor space is given by 

M = ^ [{as + bs) + {as - bs)(Ti ■ h(T2 ■ n + (a^ - hs)(Ti ■ 7x10-2 ■ rh 

+ {as - bs)(Ti ■ i(T2 ■ i + es{ai + 0-2) ■ n (24) 

where 

Z = ^-t^, rn = ^-^, n=^^ (25) 
\p' + p\ \p' — P\ \P ^ P'\ 

for p and p' as the initial and final momenta of the proton. The first term is the central 
contribution, the next three terms are double-spin-fiip terms and the final term is a single- 
spin-flip term. We can determine the sensitivity of the (e, e'p) observables to these terms by 
determining the Saclay amplitudes a^, bg, Cg, dg and from the helicity matrix elements as 
described in the Appendix, setting some of the amplitudes to zero and then transforming 
the result to give new invariant amplitudes for the Fermi form. A common approximation 
is to use only the central part of the amplitude generated from a prescription involving the 
total cross section. 

The contribution to the current matrix element given by (I2T1) . represented by Fig. [2)d, 
involves a principal value integral over off-mass- shell momenta for one leg of the scattering 
amplitude. The proton propagator for this leg contains only the positive energy contribution. 
Determination of the off shell behavior of the scattering amplitude requires a dynamical 
model of the amplitude. Such a model is not currently available to us in the range of 
energies required for the experiments being performed at Jefferson Lab. In order to estimate 
the possible effects of this contribution to the current matrix elements, we use a simple 
prescription for the off-shell behavior of the amplitude. Although additional invariants are 
possible when the nucleon is allowed to go off shell, we keep only the forms in (!23|l . The 
center-of-momentum angle is calculated using 
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FIG. 3: Off-shell nucleoli form factor for Aat = 1.0, 1.1 and 1.2 GeV. 



cos e = (26) 

The invariants are then replaced by 

Ti{s, t) Ti{s, t, u)Fn{s + t + u- 3m^) (27) 

where 

2 (A^ - 

(p^ — m^)^ + (A^ — m^Y 

and the J-'i{s,t,u) are obtained from interpolation of the on-shell invariant functions with 
the center-of-momentum angle obtained from (!26|) . The form factor (l28l) was used as a cutoff 
in calculating the Gross vertex function used in this paper with A^v = 1.675 GeV. However, 
there should be an intrinsic fall off of the scattering amplitude due to the dynamics of the 
scattering which would be expected to be faster than provided by this cutoff mass. Figure 
[3] shows the off-shell form factor for cutoff masses of Aat = 1.0, 1.1 and 1.2 GeV. In the 
absence of a dynamical model of the scattering amplitudes, the effect of possible off-shell 
contributions on various observables can be reasonably estimated by using cutoff masses in 
this range. 

The contribution to the current matrix elements from fl22]) . represented by Fig. con- 
tains the effect of negative energy propagation of the off-shell leg of the scattering amplitude. 
Since the denominator of this part of the propagator will be large compared with that of 
the positive energy part at large momentum transfers, we neglect this contribution for the 



present. This is the only approximation involved in our calculation. Note, however, that by 
dropping this contribution the current matrix elements are no longer covariant. We have 
chosen to calculate the matrix elements in the laboratory frame. 



Differential Cross Section 



The D(e, e'p) cross section for unpolarized deuterons and protons in the lab frame can 



be written as 



34 



35| 



de'dflgdQp y ^ 



C" Mott free 



TTlp TTln Pp 

{vlRl + VtRt + VttRtT cos 2(j)p + VltRlT cos (j)p) 
+hVLT'RLT' sin (j)p , 



(29) 



where M^, mp and are the masses of the deuteron, proton and neutron, pp = pi and Qp 
are the momentum and solid angle of the ejected proton, e' is the energy of the detected 
electron and fig is its solid angle. The helicity of the electron is denoted by h. The Mott 
cross section is 

(7 Mott = 



and the recoil factor is given by 

free 

The kinematic coefficients vk are 



1 + 



acos{9e/2) 
2e sin2(0,/2) 

upp — Epq cos Op 



MdPp 



= -T 



4 



r 

, , 2 

vr = ^ + tan - 



vtt = — 



2g2 



Oe 



vlt' = 1= — tan — 

V2g2 2 



(30) 



(31) 



(32) 
(33) 
(34) 

(35) 
(36) 
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If the response tensor is defined as 

W^"^^ {PlSr,P2S2\J^\P\dY {PlSl,P2S2\r\P\d) (37) 

the response functions Rk are defined by 

Rt = W^^ + W'^'^ = wi,! + w.i _i 

RTTCOs2(f)p = W^"^ - W^^ = 2^{wi^_i) 

Rlt cos (l)p = 2V2U{W^^) = -23?(w? - w°_i) 

Rlt' sin(/)p = -2V2Q(W'^^) = -23?(w° + -u;"!) , (38) 

WX'^,X^^^ {PlSi;P2S2\J\'^\PXd)*{PlSi;P2S2\Jx^\P^d) (39) 

I] (Pi5i;P2S2| J°|PAd)*(piSi;p2S2| Ja, li'Ad) (40) 



where 



and 



.T, 3 



with 



J±i = T^{J' ± J') (41) 
For our calculations, we have chosen the following kinematic conditions: the z-axis is 
parallel to q, and the missing momentum is defined as = Q — Pi — P2 — Pn- 



B. Asymmetries 

The representation of the cross section in terms of response functions is due to the mixed 
polarization of the virtual photon which varies with the electron kinematics and polarization. 
The transverse-transverse response function is the result of interference between the = 
±1 helicity states while the longitudinal-transverse response functions Rlt and Rlt' are 
the result of interference between the deuteron charge and two linear combinations of the 
A^ — ±1 helicity states. As an alternative to a complete separation of the cross section 
into response functions, the interference response functions can be accessed through linear 
combinations of differential cross sections to produce three interference asymmetries defined 
as 

Att = ""^^^^^ (42) 
vlRl + vtRt ^ ' 
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^ ^ ao(0°)-ao(180°) ^ vltRlt .^3. 

ao(0°) + (To(180°) vlRl + vtRt + vttRtt 

and 

^ ^ a+i(90°)-a_i(90°) ^ ^LT'i^LT' 

a+i(90°) +(T_i(90°) vlRl + vtRt -VttRtt' 

where, for conciseness, 

Note that while Alt can be obtained by measuring protons in the electron scattering plane 
symmetrically about the direction of the three-momentum transfer, the asymmetries Att 
and Alt' require measurements to be made out of the scattering plane. The asymmetry Alt' 
is defined as an electron single spin asymmetry and can, therefore, be easily obtained by 
flipping the beam helicity. While Rl and Rt are independent of photon-helicity-dependent 
phases, the interference response functions are not. As a result, the interference response 
functions can be very sensitive to phase differences generated by non-nucleonic currents and 
final state interactions. This is particularly true of Rlt' which can be shown to be zero in the 
PWIA. The interference response function, Rlt-, is very sensitive to the relativity included 
in the current operator, due to the various interference contributions from the charge and 



transverse current operators 



36|. 



The observable Alt' has recently been measured in Jefferson Lab's Hall B 



17|. Due to 



the large solid angle coverage in Hall B, and the averaging over 0p, the transverse-transverse 
interference response that is multiplied with a factor of cos(2</)p) in the cross section, drops 
out of the measured asymmetry Alt' js^]- Therefore, from now on in this paper, we calculate 

^HallB^ VLT'RlT' 

vlRl + vtRt ^ ' 

The difference between the asymmetry calculated with and without Rtt is very small in 
practice, due to the small size of the transverse-transverse response. 



III. RESULTS 

In this section, we discuss our numerical results for several observables. We will investigate 
the effect of the final state interactions (FSIs), and in particular, we will point out the 
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contributions of spin-dependent FSIs, both the single-spin-flip contributions and the double- 
spin-flip contributions. We will also discuss the relative importance of on-shell and off-shell 
contributions to FSI. 



A. Differential Cross Sections 



1. Momentum Distributions 



The most natural observable to investigate is the differential cross section. In Fig. [Hour 
calculations are compared to the published Die, e'p) data from Jefferson Lab's Hall A. The 

n 

data from Ulmer et al. |2(| are shown together with our PWIA, on-shell FSI, and full FSI 
curves as a function of the missing momentum p^- For all of the calculations in the paper 
we use the MMD 38| proton electromagnetic form factors. The data are presented as a 



reduced cross section, which is deflned as 



^reduced — , j-, • ) 

dileidilj^atjf.' aepTnpTnnPp 

At lower missing momenta, the effects of FSIs are small. For missing momenta larger 
than pm ~ 110 MeV, the FSI effects become visible. First, the PWIA curve is slightly above 
the FSI results, then, at pm ~ 300 MeV, the FSI curves become larger than the PWIA 
contribution. The agreement with the data is quite nice overall. The off-shell FSI is small, 
and leads to flnal results a little below the data at larger missing momenta. This is a sensible 
result, as we do expect meson exchange currents (MECs) to play a role at the relatively low 
at which these data were taken. Indeed, in [2Q], the calculation by Arenhoevel |39| 
agreed with the data at large missing momenta after MEC contributions were included; the 
FSI-only calculation was a bit below the data. 

We also have added a panel with a linear plot of just the low missing momentum data. 
The full FSI curve and the on-shell FSI curve coincide in this region. It has been observed in 
several previous D{e, e'p) measurements that at very low missing momenta, the calculations 
are somewhat above the data. This is quite puzzling as at these low missing momenta, 
effects like FSIs, MECS etc are supposed to be small and well under control. For a nice 



compilation on this topic, see [2l|]. In [201], Fig. 1 shows the deviation of the reduced cross 
section data and the calculations. Here, we observe the same type of deviation at very low 
Pm < 0.35 MeV. The largest discrepancy appears at Pm = 15 MeV, where our calculation 
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FIG. 4: The reduced cross section for a beam energy of 3.1095 GeV, = 0.665 GeV^, xbj = 0.964, 



and (pp = 180°. The data are from 



2C| 



overpredicts the data by 15%. Overall, comparing with the previous results, our low missing 
momentum results seem to be an improvement, even though the discrepancy has not been 
fully removed. One main difference between the calculation presented in this article and the 
calculations in 20] is the fully relativistic approach we take here. 



In Figure O the upper panel shows the cross section as a function of the missing momen- 
tum, with a beam energy of 5.5 GeV, = 2 GeV^, xbj = 1, and (pp = 180°. The choice 
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FIG. 5: Top panel: the differential cross section for a beam energy of 5.5GeV, = 2 GeV^, 
XBj = 1, and 0p = 180° is shown in PWIA (dotted curve), with on-shell FSI (solid curve), and 
with on-shell and off-shell FSIs (dashed curve), as a function of the missing momentum. Bottom 
panel: the ratio of the off-shell calculations with varying cut-off Ajv = 1 GeV (short-dashed), 
1.1 GeV (dash-dotted), and 1.2 GeV (long-dashed) to the on-shell FSI calculation, in the same 
kinematics as the top panel. 
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of Bjorken-x implies roughly quasi-free kinematics. The azimuthal angle of the detected 
proton has been chosen to maximize the cross section. The PWIA and on-shell FSI curves 
are almost identical at very low missing momenta, up to 0.05 GeV. Then, the FSI curve 
reduces the cross section compared to the PWIA result, roughly from = 0.05 GeV to 
Pm = 0.33 GeV. For larger missing momenta, there is a marked increase in the differential 



cross section when FS 
other calculations [23, 



is included. These results are quite typical and have been seen in 
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26l . |40| . The differential cross section decreases by several orders 
of magnitude with increasing missing momentum, and a small reduction due to FSI at lower 
missing momenta and larger differential cross section can lead to a very large increase at 
larger missing momenta and smaller differential cross sections. The inclusion of the off- 
shell FSI contributions leads to a slight reduction of the cross section for medium missing 
momenta. 

The lower panel shows the ratios of the off-shell FSI calculations with A^r = 1.0, 1.1 
and 1.2 GeV to the on-shell FSI. The off-shell effects are small for small pm, but become 
increasingly large as pm increases. These effects are not particularly sensitive to the cutoffs 
chosen here, but must be quite sensitive to lower cutoff masses. 

In Fig. [6l we investigate the role played by the various parts of the proton-neutron 
scattering amplitude contributing to (!23|) . The dashed line shows the results without the 
three double-spin-flip terms in the amplitude. One can see that for missing momenta from 
0.2 GeV to 0.45 GeV, the double-spin-flip contribution is quite relevant. Its omission in 
this region leads to a significantly smaller cross section. The single-spin-fiip - or spin-orbit 
- part of the proton-neutron scattering amplitude becomes relevant only at higher missing 
momenta than the double-spin-fiip terms. From roughly p^ = 0.3 GeV on, omitting the 
spin-orbit contribution leads to a decrease in the differential cross section. 

This clearly shows that while it is possible to parameterize the central part of the A^A^ 
scattering amplitude, and to reproduce the A^A^ cross section data this way, this type of 
parametrization effectively includes some physics that stems from the spin-dependent parts 
of the A^A^ amplitude. Here, it is interesting to see the influence of the spin-dependent parts 
of the A^A^ amplitude on the unpolarized cross section. While the logarithmic scale necessary 
for the momentum distribution conceals the effects, it is important to note that the relative 
importance of the spin-dependent FSI contributions changes with missing momentum. We 
will return to the effects of spin-dependent FSI with Fig. [10 
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FIG. 6: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, xbj = 1, and 
<j)p = 180° is shown calculated with on-shell FSI, as a function of the missing momentum. The 
solid line shows the result calculated with the full A^A^ scattering amplitude, the dashed line shows 
the result without the double-spin- flip terms of the A^A^ scattering amplitudes, and the dotted line 
shows the result with the central A^A^ amplitude only. 

In Fig. U\ we display our results for the same four- momentum transfer, but higher xbj = 
1.3. At this value of xr,-, and larger values, strong short-range pn correlations have been 

nn 

reported by inclusive electron scattering experiments on deuterium [15|, |41|]. The deviation 
of the on-shell FSI from the PWIA is small for medium missing momenta, and seems to 
disappear altogether for missing momenta between 0.5 GeV and 0.6 GeV. However, the 
off-shell contribution to the FSI gains in relevance for larger missing momenta, and leads to 
a significant increase over the PWIA results at pm > 0.4 GeV. The lower panel again shows 
the ratio of off-shell FSI to the on-shell result for three values of the off-shell cutoff. The 
importance of the off-shell FSI here is larger than for xbj = 1 (as discussed above when 
considering Fig. [5] and the sensitivity to the value of the cutoff is much greater. This is the 
expected behavior, as the deviation from xbj = ^ corresponds to a deviation from the quasi- 
elastic kinematics, and stresses the off-shell region more. A recent new proposal 18| suggests 
a measurement of the cross section at somewhat larger = 3.5 GeV^, but the same value 
of XBj = 1-3 and a beam energy of Eheam = 5.25 GeV. The results of our calculation for 
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FIG. 7: Top panel: the differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, 



XBj 



1.3, and 



180° is shown in PWIA (dotted curve), with on-shell FSI (solid curve), and 



with on-shell and off-shell FSIs (dashed curve), as a function of the missing momentum. Bottom 
panel: the ratio of the off-shell calculations with varying cut-off Ajv = 1 GeV (short-dashed), 
1.1 GeV (dash-dotted), and 1.2 GeV (long-dashed) to the on-shell FSI calculation, in the same 
kinematics as the top panel. 



18 



these kinematics are similar to what is displayed above in Fig. [71 The on-shell FSI in this 
case deviates a bit more from the PWIA result than for the kinematics displayed here. The 
off-shell contribution is just as significant in the proposed kinematics. 

2. Angular Distributions 

Bjorken x 




FIG. 8: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, pm = 0.4 
GeV, and (pp = 180° is shown in PWIA (dotted curve), with on-shell FSI (solid curve) FSI (dashed 
curve), as a function of the angle of the missing momentum. The top axis shows the corresponding 
values of Bjorken-x. 

In Figure [8], we show the cross section as a function of the angle 6m of the missing 
momentum, with fixed = 2 GeV^ and Pm = 0.4 GeV. The beam energy and azimuthal 
angle of the proton are the same as for the momentum distribution graphs. The angular 
distribution of the cross section shows much less variation in the magnitude, and therefore 
can be shown on a linear plot, allowing for a better look at the relevance of various parts of the 
cross section. While the region beyond 6m = 80° is kinematically accessible to experiment, 
a calculation in this region requires the knowledge of the NN scattering amplitude above 
lab kinetic energies of 1.3 GeV in the pn system, which are not available from SAID. We 
therefore stay below 80°. As is obvious from the plot, the most interesting features of the 
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calculation are located below this angle: while the PWIA results are gently sloping upwards 
and then downwards for angles larger than 20°, the FSI results initially follow this behavior, 
but then show a pronounced peak at around 6m = 70°. This value corresponds to Xbj = 1, 
i.e. quasi-free kinematics for the knocked-out nucleon. For the lower angles, the FSI simply 
leads to a reduction in the cross section, but the shape is unchanged. For larger angles, 
around xbj = 1, the diffractive nature of the FSI leads to a redistribution of strength from 
smaller missing momenta, causing a large peak. 



If we consider the ratio of FSI to PW 
paring various methods of calculation 21 



A cross section, as is sometimes done when com- 



23| . this ratio peaks at 70°, too. Our calculation 



clearly shows the same shift from a peak at 90°, as seen in Glauber theory calculations 



to a lower angle, as seen in the Generalized Eikonal Approximation (GEA) 



42|, 



and the 



diagrammatic approach of Laget 



26|. 
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FIG. 9: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, pm = 0.4 GeV, 
and (j)p = 180° is shown with on-shell FSI (solid curve), and with on-shell and off-shell FSIs at 
various cut-offs, as a function of the angle of the missing momentum. 

For the angular distribution, the influence of the off-shell FSI cut-off is clearly visible and 
can be studied easily. One does not expect a large contribution from far off-shell nucleons. 
The cut-off we use here serves to impose that constraint. For the kinematics displayed in 
Fig. [9], we investigate the effects of various cut-off values. The cut-off at 1 GeV leads to a 
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very small increase at small angles and a very small decrease at large angles, but the overall 
result is hardly different from the on-shell FSI only result. The primary effect for larger 
cutoff masses is to fill in the minimum in the on-shell result from 30° to 65°. 
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FIG. 10: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, pm = 0.4 
GeV, and (pp = 180° is shown calculated with on-shell FSI, as a function of the angle of the missing 
momentum. The solid line shows the result calculated with the full A^A^ scattering amplitude, the 
dashed line shows the result without the double-spin- flip terms of the A^A^ scattering amplitudes, 
and the dotted line shows the result with the central A^A^ amplitude only. 

In Fig. [TOl we show the effects of the various spin-dependent parts of the pn scattering 
amplitude. It is very interesting to observe that in the peak region, the contribution of 
the spin-dependent FSIs (both single and double-spin-flip) amounts to about one quarter 
of the cross section. This is certainly a rather significant contribution. The contribution 
from the single-spin-flip term and the three double-spin-flip terms is about equal in the 
peak region. The figure shows that the double-spin-flip contribution to the cross section 
at small angles is almost negligible. It becomes noticeable at 9^ ~ 40°, and then leads to 
a sizable increase of the differential cross section in the peak region. The omission of the 
single-spin-flip contribution leads to a noticeable reduction in the cross section for all angles. 
The effect is most pronounced in the peak region and in the very shallow dip just before the 
peak region. 
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FIG. 11: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, pm = 0.4 GeV, 



and 



180° is shown with on-shell FSI as a function of the angle of the missing momentum. 



The solid curve shows the contribution of the full A^A^ scattering amplitude, the other curves show 
the results for just one invariant term in the A^A^ scattering amplitude. 



40|. Even 



To this point, we have considered the different contributions of the spin-dependent parts 
of the amplitude, i.e. of the A^A^ amplitude split up following the Saclay convention ( 12^ . 
Using the Saclay formalism with its classification according to the spin-dependence is quite 
useful, as it allows one to understand the new contributions from different parts of the current 
operator when adding the single-spin-fiip term and the double-spin-fiip terms, see 
though we did not rewrite our current operator to distinguish e.g. between magnetization 
current and convection current, seeing the A^A^ scattering amplitude in terms of its spin- 
dependence is a very natural view point, and allows for a certain intuitive understanding of 
the numerical results. 

It is also interesting to investigate the A^A^ amplitude in terms of the five Fermi invariants, 
which are so practical for actual calculations. /^From (]A8|A22|) . it is clear that every invariant 
amplitude contains several different, spin-dependent pieces. In Fig. [TTl we show the results of 
the calculation with on-shell FSI if only one of the five invariants is included. For comparison, 
the solid lines depicts the result obtained with the full A^A^ amplitude. The result obtained 
with just the pseudoscalar part is slightly above the full result for smaller angles, and 
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continues smooth and almost straight towards larger angles. It does not exhibit a peak 
structure at large angles. The tensor and axial vector contributions are fairly close to the 
pseudoscalar contribution at large angles, and neither exhibits a peak at large angles. At 
small angles, however, these two contributions lead to a new peak, much larger than the 
original peak at large angles in the full result. The scalar and vector contributions show 
peaks at small angles, and a large peak at large angles. These results already show that there 
are large interference effects present between the various invariant amplitudes. There is no 
straightforward and intuitive explanation available for why these contributions interfere in 
such a way. 
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FIG. 12: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV'^, pm = 0.2 GeV, 
and (pp = 180° is shown in PWIA (dotted curve), with on-shell FSI (solid curve), and with on-shell 
and off-shell FSIs (dashed curve), as a function of the angle of the missing momentum. 

In Fig. [121 show the angular distribution of the differential cross section for a lower 
missing momentum, = 0.2GeV. The lower missing momentum implies that the limiting 
value of 1.3 GeV lab energy for the A^A^ system is reached at larger angles than for pm = 
0.4 GeV. While the PWIA curve here is very similar in shape to the curve at the higher 
missing momentum value, the FSI curve looks rather different now: instead of a fairly sharp 
peak around 9m = 70°, we now observe a broad, shallow dip at larger angles. At lower 
missing momenta, the FSIs lead to a reduction in the cross section. Part of this strength 
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is redistributed to larger missing momenta, as discussed above. Also including the off-shell 
FSIs has no effect at very small angles, roughly below 20°, and then tends to shift the overall 
result towards lower angles. 




FIG. 13: The differential cross section for a beam energy of 5.5 GeV, = 2 GeV^, pm = 0.2 
GeV, and (pp = 180° is shown calculated with on-shell FSI, as a function of the angle of the missing 
momentum. The sohd line shows the result calculated with the full A^A^ scattering amplitude, the 
dashed line shows the result without the double-spin- flip terms of the A^A^ scattering amplitudes, 
and the dotted line shows the result with the central A^A^ amplitude only. 

In Fig. [131 the influence of the different kinds of FSIs is shown at low missing momentum. 
Switching off the double-spin-flip contribution leads to a small reduction in the cross section 
for medium and large angles, roughly 5% in the region of the shallow dip. Switching off the 
single-spin-flip term, too, changes practically nothing. The FSIs in this kinematic region 
are overall smaller than for higher missing momenta. The influence of spin-dependent FSIs 
is smaller here, too. However, it is interesting to note that the double-spin-flip terms are 
actually more relevant here than the single-spin-flip terms. We have observed this already 
when discussing the momentum distributions shown in Fig. O 
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B. Asymmetries 



Figure [H] shows the interference asymmetry Att for = 2.0 GeV and Ebeam = 5.5 
GeV. The upper panel is for xbj = ^ and the lower for xbj = 1-3. 





FIG. 14: Top panel: the asymmetry Att is shown for xbj = 1 in PWIA (dotted hne), with on-shell 
FSI (sohd line), and with off-shell FSI for various values of the cutoff - Ajv = 1-0 GeV (short-dashed 
line), Atv = 1-1 GeV (dash - dotted line), and A^v = 1-2 GeV (long-dashed line). Bottom panel: 
the asymmetry Att is shown for xbj = 1-3. The curves have the same meaning as in the top 
panel. 

In both cases, the FSI result in a substantial change from the PWIA result, both with 
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respect to size and shape of the asymmetry. For xbj = 1, the variation with the cutoff mass 
for the three values shown here is small, but it is somewhat larger for xbj = 1-3 reducing 
the asymmetry almost to zero for A^r = 1.2 GeV around Pm = 0.5 GeV. 
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FIG. 15: Top panel: the asymmetry Alt is shown for xbj = ^ in PWIA (dotted line), with on-shell 
FSI (solid line), and with off-shell FSI for various values of the cutoff - A^v = 1-0 GeV (short-dashed 
line), Atv = 1-1 GeV (dash - dotted line), and Ajv = 1-2 GeV (long-dashed line). Bottom panel: 
the asymmetry Alt is shown for xbj = 1-3. The curves have the same meaning as in the top 
panel. 



Figure [T^ shows results for the interference asymmetry Alt for the same kinematics as in 
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the previous figure. Again the sensitivity to final state interactions is substantiaL Sensitivity 
to off-shell contributions is relatively modest at xbj = ^ but is much larger at xbj = 1-3. 

The single-spin asymmetry A^t' is shown in Fig. for the same kinematics. This 
asymmetry is identically zero in the PWIA, so the FSI are responsible for any asymmetry. 
A more detailed discussion of FSI effects is given below for kinematics relevant to a recent 
experiment. Here, in these kinematics, we focus on comparing the behavior of the three 
asymmetries, in the same kinematics. 

At XBj = 1 there is very little sensitivity to off-shell contributions, but it is large at 
xsj = 1-3. As stated above, this behavior is expected as xsj = 1 corresponds to quasi-free, 
i.e. on-shell, kinematics, while xbj = 1-3 probes nucleons that are much more off-shell. 
From the plots for xbj = 1.3, one can see clearly that the off-shell contribution to the NN 
scattering amplitude introduces a certain amount of ambiguity, especially at medium to 
high missing momenta. While great progress has been made on the experimental side, with 
measurements at very high missing momenta, there clearly are some theoretical uncertainties 
in these kinematic regions. 

We are in the fortunate situation that the asymmetry Alt' has been measured over a 
wide range of kinematics, from very low four-momentum transfers ~ 0.2 GeV^ up to 
medium values of ~ 1 GeV^ 17|. The data are currently under analysis. In this range, 
the proton-neutron scattering amplitudes from SAID are available, so there are no limits to 
our ability to calculate for these kinematics. 

Here, we discuss our results for two representative kinematics: = 0.5 GeV^, and 
Q2 = 1.0 GeV^ In both cases, we assume a beam energy of 2.558 GeV. Fig. [T71 shows our 
results for = 1 GeVl The PWIA result is zero, and not shown on the plot. The full, 
on-shell FSI result starts out negatively, dips around pm ~ 0.25 GeV, and then increases and 
changes sign around pm ~ 0.38 GeV. Then, the asymmetry peaks around pm ~ 0.5 GeV, 
and then decreases and changes sign again. The second zero of the asymmetry occurs at 
Pm ~ 0.7 GeV. One can see that the double-spin-flip contributions to the FSI are not very 
relevant: they just lead to small modiflcations in the dip and peak regions. The spin-orbit 
(i.e. the single-spin-flip) contribution is very important, though. Switching it off so that 
only the central FSI remains leads to a completely different picture: the asymmetry is tiny, 
and remains positive for large missing momenta. 

In Fig. [TSl we show the corresponding results for the spin-dependence of the FSIs for 
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FIG. 16: Top panel: the asymmetry Aix' is shown for xbj = 1 in PWIA (dotted line), with 
on-sheh FSI (solid line), and with off-shell FSI for various values of the cutoff - Ajv = 1-0 GeV 
(short-dashed line), Ajv = 1-1 GeV (dash - dotted line), and Ajv = 1-2 GeV (long-dashed line). 
Bottom panel: the asymmetry A^^t' is shown for xsj = 1-3. The curves have the same meaning as 
in the top panel. 

lower four-momentum transfer, = 0.5 GeV^. At these kinematics, we do expect the 
influence of meson exchange currents and isobar states to be relevant. These effects are not 
included in the present calculation. However, the FSIs are crucial for Alt', and we can 
investigate them within our model. The full calculations are qualitatively very similar at 
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FIG. 17: The LT' asymmetry for a beam energy of 2.558 GeV and four-momentum transfer 
= 1 GeV^ is shown calculated with on-shell FSI, as a function of the missing momentum. The 
solid line shows the result calculated with the full A^A^ scattering amplitude, the dashed line shows 
the result without the double-spin-flip terms of the A^A^ scattering amplitudes, and the dotted line 
shows the result with the central A^A^ amplitude only. 

both values, showing a negative dip around Pm ~ 0.25 GeV and then an increase into 
positive values, with a peak around pm ~ 0.5 GeV. A very interesting difference, however, is 
the size of the contribution of the double-spin-flip terms to the FSI. While their influence is 
small, almost negligible, at the higher value, it is quite significant for the low value: 
the double-spin-flip terms serve to partially fill in the negative dip, and are also responsible 
for pushing the asymmetry back towards positive values. 

C. FSI Details 

One obvious difference between the calculation presented here and the traditional Glauber 
and generalized eikonal approximation (GEA) is the evaluation of the argument of the 
nucleon-nucleon scattering amplitude for the FSIs. As described in section[TTl we evaluate the 
five terms of the NN scattering amplitude fl23l) at the values of the Mandelstam variables s and 
t computed from the particular kinematics. In Glauber and GEA settings, one typically finds 
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FIG. 18: The LT' asymmetry for a beam energy of 2.558 GeV and four-momentum transfer 
= 0.5 GeV^ is shown calculated with on-shell FSI, as a function of the missing momentum. 
The solid line shows the result calculated with the full NN scattering amplitude, the dashed line 
shows the result without the double-spin-flip terms of the A^A^ scattering amplitudes, and the 
dotted line shows the result with the central NN amplitude only. 

expressions where the NN scattering amplitude is evaluated assuming a purely transverse 
momentum transfer when evaluating t, even though the longitudinal momentum transfer 
is taken into account in the GEA. This is typically denoted with expressions like fNN{k±)- 
With the kinematics variables as defined in Fig. [T|, the Mandelstam t is given by t = 
(Mrf — + uj — Ep-^)^ — (pi + k2 — q)^, whereas assuming a purely transverse momentum 
transfer implies: t± = — (pi,± + k2,±)^. Using this we can define the cm scattering angle as 



In Fig. [19], we show the ratio of the transverse-momentum approximation to the full, 
on-shell calculation. The kinematics are identical to the kinematics used for Fig. (5) a beam 
energy of 5.5 GeV, = 2 GeV^, xbj = 1, and (pp = 180°. Up to missing momenta of 0.4 
GeV, the two approximation works well, leading to small deviations of less than 5%. Beyond 
0.4 GeV, the deviation from the full result grows, and for missing momenta larger than 0.6 
GeV, the quality of the approximation deteriorates quickly. 



cos 6* = 1 + 



2t± 



(48) 



s — 4m^ 
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We have performed the same calculation for the angular distribution shown in Fig. [81 
Here, we found that the approximation does well - the deviations are less than 5% for any 
angle. This corresponds to our findings for the momentum distribution. 

As an illustrative example for the other observables discussed in this article, we also show 
our results for the asymmetry A^t' in Fig- 1201 The kinematics correspond to Fig. [T71 we 
have used a beam energy of 2.558 GeV and four-momentum transfer = 1 GeV^. The 
spike seen in the ratio around pm = 0.4 GeV stems from the sign change in the asymmetry, 
they are not relevant to our discussion. Again, we see that the approximation is doing 
well up to missing momenta of roughly 0.4 GeV. Then, the approximation considerably 
overestimate the full result, and for very large missing momenta, > 0.7 GeV, it even fails 
to reproduce the correct sign of the asymmetry. It is interesting to see that the effects of 
the approximation are visible even in a quantity that is a ratio of quantities that are both 
affected by the approximation. 

In conclusion, approximating the argument of the NN scattering amplitude with the 
popular "transverse momentum transfer only" works well up to missing momenta of 0.4 
GeV for various observables. For missing momenta higher than that, the approximation 
becomes questionable. 



IV. SUMMARY AND OUTLOOK 



In this paper, we have presented a fully relativistic D{e,e'p) calculation in impulse ap- 
proximation. We have steered clear from a number of common, simplifying assumptions. 
The only approximation made in this paper is neglecting the negative-energy contributions 
to the propagator of the off-shell nucleon. These contributions can realistically expected to 
be very small compared to the positive-energy contributions. 

We have used a parametrization of experimental NN data from SAID to describe the full 
pn scattering amplitude for the final state interaction. This leads to certain limits in the 
kinematics we can access, as these parametrizations are available only for lab kinetic energies 
of 1.3 GeV or less. In our calculations, we have investigated the effects of the different 
contributions to the NN scattering amplitude: the central, spin-orbit, and double-spin-flip 
parts, using the Saclay formalism to describe the different contributions. Many other fine 



calculations using the generalized eikonal approximation 23|, UM, 1251] or a diagrammatic 
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FIG. 19: The ratio of the differential cross section evaluated with the full A^A^ argument, t, and 
with the purely transverse momentum, t±. The kinematics are identical to Fig. [5j beam energy 
of 5.5 GeV, = 2 GeV^, xbj = 1, and 0p = 180°. The calculation was performed using on-shell 
FSI. 



approach 26|] use the central part only. While the central part of the amplitude is clearly 
dominant in almost all observables, the spin-orbit and double-spin-fiip parts do contribute 
visibly to the cross section in the peak area of the angular distribution, increasing the peak 
height roughly by a quarter. For the out-of-plane asymmetry Alt/, which is non-zero only 
in the presence of FSIs, the spin-orbit part is clearly the most relevant. Depending on the 
kinematics, the double-spin-flip can also play a relevant role for this asymmetry. 

We also showed the different contributions of the NN amplitude in terms of the five 
invariant amplitudes. Interestingly, they are all relevant, and a lot of interference effects 
contribute to the full result. It is not possible to identify a single, dominant contribution in 
this framework for the description of the NN amplitudes. 

In the spirit of avoiding all unnecessary approximations, we have used the full argument 
for the calculation of the NN scattering amplitudes. In Glauber theory and its variants, one 
often encounters the assumption of a transverse momentum transfer only, and this changes 
the value of Mandelstam t. We have investigated the validity of this assumption, and found 
that it is a very good approximation for missing momenta up to 0.4 GeV. For higher missing 
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FIG. 20: The ratio of the asymmetry Ait' evaluated with the fuh A^A^ argument, t, and with the 
purely transverse momentum, t±. The kinematics are identical to Fig. [171 The calculations were 
performed using on-shell FSI. 

momenta, the quality of this approximation deteriorates quickly, and it should probably not 
be used. 

We have also compared the influence of the off-shell FSI contributions to the on-shell 
FSI contributions. The former are expected to not be too large, and they require some 
interpolation of on-shell amplitudes and the introduction of a regulator function to suppress 
very far off-shell contributions. The off-shell FSI contributions tend to be negligible to small 
for lower missing momenta, Pm < 0.4 GeV for any observable. Beyond that, their importance 
varies depending on the specific kinematics: the off-shell FSI is very small for the momentum 
distribution calculated for the quasi- free value xbj = 1, but it is large for xbj = 1-3. This 
pattern was observed both for the cross section and the asymmetries Att, ^lt? and A^j". 
The size of the off-shell contribution does depend on the chosen cut-off, with a larger cut-off 
admitting a sometimes much-larger contribution. The main purpose of showing figures with 
the ratios of off-shell calculations with different cut-offs to the on-shell result is to identify 
"safe" kinematics and observables, where the off-shell FSI contributions are definitely small. 
In regions where they are relevant, a certain amount of theoretical uncertainty cannot be 
avoided, until reliable and believable models of the off-shell NN interaction at the relevant 
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energies are developed. This is particularly relevant for the interpretation of new data taken 
at high missing momentum at Jefferson Lab. 

The current calculations will be applied or already have been applied to the forthcoming 
experimental data from Jefferson Lab 1^, , and calculations for the BLAST data from 



MIT Bates 



16| are planned. Our calculation also does a nice job of improving the agreement 



with the low missin g rn omentum data of Ulmer et al. 20|, even though the "low missing 



momentum puzzle" 



2l| is not completely resolved. 



Logical next steps for enhancing our calculations are the inclusion of meson exchange 
currents, and isobar states. 
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APPENDIX A: REPRESENTATIONS OF THE NN AMPLITUDES 

The invariant functions J^i{s,t) can be obtained from scattering data. For example, the 
helicity amplitudes are defined as 

A<A;A^;AiA2 = (MA;(p'l))a(MA'2(P2))b^ab,crf(MAi(Pl))c(^^A2(P2))d (Al) 

where ux{p) is the helicity spinor for helicity A. The helicity matrix elements for scattering 
in the center-of-momentum frame can be obtained from the program SAID for laboratory 
kinetic energies of up to 1.3 GeV. If the amplitudes are extracted in units of fm, the helicity 
amplitudes consistent with the conventions used here are related to the SAID amplitudes 
by 
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Parity, time-reversal and particle interchange symmetries can be used to show that there 
are only five independent helicity amphtudes defined as 



a 
h 
c 
d 
e 



A^l,l;l,l 

A^l,l;-1,-1 
A^l,-l;-l,l 



(A3) 
(A4) 
(A5) 
(A6) 
(A7) 



Using (1251) in flAip and solving for the invariant functions gives 
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where 
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s sin 6* 
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fs sin 6 

m^[2m^{l + cos 6*) - s(3 + cos 6*) 
s(l + cos 6*) 
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s(l + cos6') 

The Saclay amphtudes are given in terms of the hehcity amphtudes by 
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(A9) 
(AlO) 
(All) 

(A12) 

(A13) 

(AM) 

(A15) 
(A16) 

(A17) 

(A18) 

(A19) 

(A20) 

(A21) 



(A22) 



Equations flASP and flA22l) can then be used to obtain the transformation from the Saclay 
amphtudes to the Fermi invariant functions. 
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